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SECTION  1 


INTRODUCTION 


An  analysis  was  developed  in  this  Northrop/ AFWAL  program 
to  predict  the  strength  of  bolted  composite  structures.  This  report 
presents  details  of  the  developed  analysis,  sample  predictions,  and 
a  discussion  on  its  validity  and  its  application  to  structural 
design. 


.  Prior  to  the  initiation  of  this  program,  the  strength  of  a 
bolted  laminate  was  analytically  predicted  using  approximate 
analyses  and  experimental  results.  The  distribution  of  the  applied 
load  among  the  fasteners  was  initially  obtained,  and  the  most 
critical  fastener  location  was  subsequently  analyzed  to  predict  the 
joint  strength.  The  fastener  load  distribution  analysis  was 
essentially  one-dimensional,  assuming  that  all  the  fasteners  in  a 
row  (perpendicular  to  the  load  direction)  carried  equal  loads.  And, 
the  load  distribution  among  the  various  rows  was  predicted  based  on 
experimentally  obtained  "joint  stiffness"  values.  The  subsequent 
strength  analysis  at  a  fastener  location  was  based  on  an  infinite 
plate  stress  analysis  and  was  incapable  of  accounting  for 
neighboring  stress  concentrators  (like  a  free  edge,  a  cut-out  or  a 
neighboring  fastener  location) . 

A  strength  analysis  was  developed  in  this  Northrop/AFWAL 
program  to  overcome  the  major  deficiencies  that  existed  at  its 
inception.  The  analysis  incorporates  special  finite  elements  into  a 
failure  analysis  procedure  that  predicts  the  fastener  load 
distribution,  the  critical  fastener  location,  the  joint  strength  and 
its  failure  mode.  Four  special  finite  elements  were  developed  using 
a  fastener  analysis  and  a  stress  analysis  that  accounts  for  finite 


laminate  planform  dimensions  (see  Reference  1) .  These  elements 
include  a  loaded  hole  element,  an  unloaded  hole  element,  a  plain 
element,  and  an  effective  fastener  element.  A  finite  element  model 
of  the  bolted  joint  computes  the  fastener  load  distribution  and 
averaged  stresses  at  each  fastener  and  cut-out  location.  The 
critical  fastener  or  cut-out  location,  the  joint  failure  load  and 
the  corresponding  failure  mode  are  predicted  based  on  these 
computations. 

The  developed  strength  analysis  has  been  programmed  to  be 
the  SAMCJ  (Strength  Analysis  of  Hultifastener  Composite  Joints) 
computer  code.  SAMCJ  requires  a  definition  of  the  geometry  and  the 
material  properties  of  the  bolted  plates  and  fasteners  as  input. 

The  presence  of  any  cut-out  is  included  in  the  finite  element  model 
as  an  "unloaded  hole"  element  linked  to  adjacent  "loaded  hole"  and 
plain  elements.  The  input  material  properties  of  the  bolted 
laminates  inolude  failure  parameters  that  are  required  by  the 
average  stress  failure  criteria.  These  are  distances  from  the 
fastener  or  cut-out  hole  boundaries,  at  selected  locations,  over 
which  stresses  are  averaged  and  compared  to  plain  laminate 
strengths,  to  predict  failure  (see  Reference  1) .  SAMCJ  computes  the 
joint  load  values  for  net  section,  bearing  and  shear-out  modes  of 
failure  at  each  fastener  and  cut-out  location.  Information 
corresponding  to  the  least  value  provides  the  joint  failure  load, 
the  critical  fastener  or  cut-out  location,  and  the  failure  mode. 

In  computing  the  fastener  load  distribution  and  the 
critical  average  stress  values  at  every  fastener/ cut-out  location, 
SAMCJ  also  accounts  for  fastener  flexibility  effects.  The  FDFA 
fastener  analysis,  developed  earlier  in  the  program  (Reference  1), 
is  used  to  compute  the  effective  fastener  stiffness,  accounting  for 
bolt  torque  and  load  eccentricity  (single  versus  double  shear 
transfer  of  the  applied  load) .  FDFA  is  employed  twice  to  compute 
the  effective  transverse  stiffnesses  of  the  fastener,  along  and 
perpendicular  to  the  load  direction.  The  effective  fastener 
stiffness  matrix  connects  the  bolted  plates  at  the  fastener 
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locations,  accounting  for  all  significant  joint  parameters. 

The  significant  improvements  offered  by  SAMCJ  over  the 
state-of-the-art  at  program  inception  are: 

(1)  SAMCJ  performs  a  one-step  analysis  that  computes  the 
fastener  load  distribution,  critical  fastener  location,  joint 
failure  load,  and  the  corresponding  failure  mode.  Hitherto, 
separate  fastener  load  distribution  and  failure  analyses  were 
performed,  requiring  a  two-step  analytical  procedure. 

(2)  SAMCJ  only  requires  the  geometric  and  material 
properties  of  the  bolted  plates  and  fasteners  as  input.  SAMCJ 
internally  computes  the  effective  transverse  fastener  stiffness 
values  that  account  for  fastener  size,  fastener  and  bolted  plate 
material  properties,  bolt  torque,  load  eccentricity  (single  versus 
double  shear  load  transfer) ,  and  the  local  three  dimensional  stress 
state  at  the  fastener  location.  Hitherto,  these  effeots  could  only 
be  accounted  for  via  experimentally  measured  "joint  stiffnesses." 
SAMCJ  eliminates  the  need  for  these  experimental  measurements,  and 
is,  therefore,  the  first  multifastener  bolted  joint  strength 
analysis  that  is  devoid  of  dependence  on  test  results. 

(3)  SAMCJ  performs  a  two-dimensional  load  distribution 
analysis,  and  predicts  the  magnitude  and  orientation  of  the  load  at 
each  fastener  location  via  components  of  the  fastener  load  along  and 
perpendicular  to  the  load  direction.  Analyses  available  at  program 
inception  only  addressed  the  row-to-row  load  variation,  or  the  axial 
components  of  the  fastener  loads,  resorting  to  a  one-dimensional 
analysis. 


(4)  SAMCJ  accounts  for  stress  concentration  interaction 
effects  that  hitherto  could  not  be  accounted  for.  This  includes  the 
effects  of  adjacent  free  edges,  cut-outs  and  proximate  fastener 
locations . 
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(5)  SAMCJ  accounts  for  tapered  bolted  plate  geometries 
that  are  commonplace  in  practical  situations. 

There  are,  however,  segments  of  the  SAMCJ  computer  code 
that  can  be  improved  beyond  their  present  capabilities.  These  are 
addressed  in  the  following  sections  of  this  report.  Nevertheless, 
the  significant  achievements  of  SAMCJ  over  the  state-of-the-art  at 
program  inception  remain  unscathed. 

SAMCJ,  by  virtue  of  the  above  qualities,  is  an  excellent 
design  tool.  It  can  be  used  to  evaluate  different  fastener  patterns 

and  to  select  among  these  for  a  specific  loading  state  at  a  bolted 
joint  location. 
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SECTION  2 


ANALYTICAL  DESCRIPTION 


This  section  presents  an  overview  of  the  strength  analysis 
of  bolted  laminates  (the  SAMCJ  computer  coda) ,  a  description  of  the 
developed  special  finite  elements,  and  the  analytical  procedure  used 
in  SAMCJ  to  predict  fastener  loads,  the  critical  fastener  or  cut-out 
location,  the  corresponding  joint  strength  and  the  failure  mode. 

The  application  of  the  developed  analysis  in  the  design  of  bolted 
laminates,  and  the  test  requirements  for  this  application,  are  also 
included, 

2.1  Overview  of  the  Strength  Analysis  of  Bolted  Laminates 

(SAMCJ) 

Ths  development  of  a  reliable  strength  analysis  is  crucial 
to  the  design  of  highly  loaded  bolted  joints  in  composite 
structures.  As  shown  in  Figure  1,  structural  loads  translate  into 
inplane  stress  resultants  (Nx,  Ny  and  Nxy)  that  transfer  from  one 
component  to  another  (skin  to  substructure,  for  example)  through 
many  fasteners.  The  SAMCJ  (Strength  Analysis  of  Multi fastener 
Composite  Joints)  computer  code  was  developed  to  analyze  such  a  load 
transfer  situation,  to  compute  the  failure  value  of  the  applied 
load,  and  to  predict  the  critical  fastener  location  and  the  joint 
failure  mode. 

A  flow  chart  of  SAMCJ  operations  is  presented  in  Figure  2. 
As  input,  SAMCJ  requires  the  user  to  specify  how  the  bolted  plates 
are  divided  into  plain  elements  and  elements  with  loaded  or  unloaded 
holes.  The  bolted  plates  are  currently  assumed  by  SAMCJ  to  be 
subjected  to  uniaxial  tensile  or  compressive  loading,  in  a  single  or 
double  shear  configuration.  The  uniaxial  restriction  on  the  load 
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Figure  1.  Schematic  Breakdown  of  the  Strength  Analysis  of  Bolted 
Joints , 


Figure  2.  Flow  Chart  of  SAMCJ  Operations. 


KKSHSKKS^ 


was  incorporated  into  the  code  because  all  the  test  specimens  in  the 
experimental  nart  of  the  program  were  only  subjected  to  uniaxial 
loads  (see  Reference  2).  This  restriction,  however,  can  be  easily 
removed  from  the  SAMCJ  code  by  modifying  it  to  require  the  user  to 
specify  the  general  biaxial  loading.  Additional  input  requirements 
for  the  SAMCJ  code  include  the  material  properties  of  the  bolted 
plates  and  fasteners,  and  the  fastener  size,  location  and  torque. 

The  material  properties  of  the  bolted  laminates  include  the  tensile 
and  compressive  failure  strains  in  the  fiber  direction  of  the 
lamina,  and  the  characteristic  distances  over  which  stresses  are 
averaged  to  predict  net  section,  shear-out  and  bearing  failures  at 
the  fastener  or  cut-out  location. 

With  the  above  input,  performs  the  following 

computations.  It  initially  generates  stiffness  matrices  for  all  the 
elements,  namely,  plain  elements,  elements  with  loaded  or  unloaded 
holes,  and  effective  fastener  elements.  The  individual  stiffness 
matrices  are  subsequently  assembled  to  obtain  the  global  stiffness 
matrix  for  the  bolted  joint.  A  1-kip  uniaxial  tensile  or 
compressive  joint  load  is  imposed  on  the  left  end  of  the  top  plate, 
in  accordance  with  the  input  instructions  (see  Figure  3) .  The  nodes 
at  the  right  end  of  the  bottom  plate  are  constrained  from 
translating  in  the  load  direction,  and  one  of  these  nodes  is  also 
constrained  in  the  transverse  direction  to  preclude  rigid  body 
translations,  (see  Figure  3).  The  solution  to  this  finite  element 
formulation  of  the  bolted  joint  provides  the  axial  and  transverse 
components  of  the  load  at  every  fastener  location,  corresponding  to 
a  l-kip  joint  load.  Also  computed  are  the  average  net  section, 
shear-out  and  bearing  stresses  at  every  fastener  and  cut-out 
location,  corresponding  to  a  1  kip  joint  load. 

SAMCJ  provides,  as  output,  the  failure  value  of  the 
uniaxial  joint  load,  the  critical  fastener  or  cut-out  location,  and 
the  joint  failure  mode.  These  are  obtained  as  follows.  The 
tensile,  compressive  and  shear  strengths  of  the  plain  laminates  are 
computed  based  on  the  input  tensile  and  compressive  failure  strains 
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in  the  fiber  direction  of  the  lamina.  The  ratios  of  the  averaged 
stresses  to  the  corresponding  plain  laminate  strengths,  at  selected 
locations  around  each  fastener  and  cut-out  boundary,  are  compared  to 
predict  the  failure  mode,  the  critical  fastener  or  cut-out  location 
and  the  joint  failure  load.  SAMCJ  predicts  net  section,  shear-out 
and  bearing  modes  of  failure  at  the  laminate  level.  In  Reference 
1,  similar  failure  predictions  for  single-fastener  joints  in 
composites  were  made  at  the  lamina  level  using  the  SASCJ  computer 
code.  Consequently,  the  failure  parameters  (characteristic 
distances  for  the  three  failure  modes)  used  with  SAMCJ  are  different 
from  those  used  with  SASCJ. 

The  incorporation  of  the  transverse  effective  fastener 
stiffness  values  provides  SAMCJ  the  capability  to  aacount  for 
fastener  flexibility,  torque,  and  load  eccentricity  (single  versus 
double  shear  load  transfer) .  The  FDFA  code,  developed  in  this 
program  (Reference  1) ,  is  used  to  compute  the  effective  faetener 
transverse  stiffnessos,  along  and  perpendicular  to  the  load 
direction.  The  effect  of  the  laminate  stacking  sequence  is  also 
accounted  for  in  this  analysis.  SAMCJ  executes  FDFA  twioe  to 
account  for  the  layup  variation  (by  90  degrees)  from  the  loading 
direction  to  the  perpendicular  direction. 

SAMCJ  accounts  for  stress  concentration  interaction 
effects  introduced  by  neighboring  cut-outs,  free  edges  and  proximate 
fastener  locations.  This  is  made  possible  by  the  use  of  the  FIGEOM 
stress  analysis,  developed  in  this  program  (Reference  1),  to 
generate  element  stiffness  matrices.  FIGEOM  accounts  for  finite 
planforra  plate  dimensions  through  a  boundary  collocation  solution 
procedure. 


SAMCJ  computes  the  magnitude  and  the  orientation  of  the 
load  at  each  fastener  location.  It  is  a  two-dimensional  load 
distribution  analysis  that  does  not  rely  on  an  experimental 
measurement  of  "joint  stiffness"  required  by  other  analyses  prior  to 
the  initiation  of  this  program.  In  a  design  situation,  many 


fastener  arrangements  can  be  analytically  and  economically  evaluated 
by  SAMCJ  to  arrive  at  the  beet  fastener  pattern  for  the  assumed 
loading  conditions. 

When  the  bolted  plates  are  tapered,  the  SAMCJ  user  can 
input  equivalent  uniform  thickness  elements  to  approximate  the 
tapering  effect.  Adjacent  elements  in  the  tapered  plate  will  have 
different  thickness  values.  This  feature  is  essential  in  the 
analysis  of  most  practical  joints. 

SAMCJ  has  been  developed  for  the  strength  prediction  of 
bolted  laminated  structural  parts.  It  currently  assumes  that  the 
selected  fasteners  preclude  fastener  failure.  Also,  it  applies  the 
same  failure  procedure  to  both  the  bolted  plates,  accounting  for  net 
section,  shear-out  and  bearing  failures  via  the  average  stress 
failure  criteria.  The  compos ite-to-metal  joints  tested  in  this 
program  (References  2  and  3)  were  designed  to  preclude  metallic 
failures.  Therefore,  validated  failure  parameters  were  not 
generated  for  the  metallic  plates.  However,  available  test  results 
in  the  open  literature  may  be  used  to  generate  these  failure 
parameters,  if  needed. 

2.2  Development  of  Special  Finite  Elements 

A  bolted  joint  region  in  a  structural  part  presents  many 
difficulties  in  performing  an  accurate  analysis.  The  local  three- 
dimensional  effects  and  the  implicit  indeterminacy  in  the  fastener 
load  distribution  add  to  the  complexity  of  the  analysis.  An 
approximate  solution  procedure  is  therefore  mandatory  in  developing 
a  strength  analysis  for  multifautener  bolted  joints.  A  finite 
element  approach  was  adopted  in  the  development  of  the  SAMCJ 
oomputer  code.  Special  finite  elements  were  developed  to  adequately 
represent  the  complex  stress  state  in  the  neighborhood  of  fasteners 
and  out-outs,  and  to  account  for  the  offective  transverse  fastener 
stiffnesses.  The  following  sub-sections  describe  the  development  of 
the  special  finite  elements. 
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2.2.1 


The  Effective  Fastener  Element 


The  effective  fastener  element  represents  the  response  of 
a  fastener  constrained  by  Its  surrounding  medium.  Its  displacement 
Is  a  function  of  many  joint  variables  which  Include  the  fastener 
size,  modulus  and  torque,  the  bolted  plate  material  properties  and 
the  layup  of  laminated  plates,  the  load  eccentricity  (single  versus 
double  shear),  and  the  geometry  of  the  joint.  In  Reference  1,  a 
finite  difference  fastener  analysis  (the  FDFA  computer  code)  was 
developed  to  account  for  the  effects  of  these  joint  parameters  on 
the  fastener  deflection.  This  was  accomplished  by  modeling  an 
Isolated  fastener  as  a  Timoshenko  beam  resting  on  an  elastic 
foundation.  A  brief  synopsis  of  the  salient  features  of  this 
analysis  Is  presented  here  for  completeness. 

Figures  4  to  6  Illustrate  the  qualitative  Influence  of 
joint  configuration  (single  versus  double  shear),  fastener  geometry 
and  properties,  and  fastener  end  constraints.  These  effects  are 
accounted  for  in  the  mathematical  representation  shown  In  Figure  7, 
and  by  the  incorporation  of  appropriate  boundary  conditions  (Figures 
8  and  9) .  The  elastic  foundation  moduli  for  the  bolted  plates  are 
plywise  uniform  in  a  laminate,  and  are  computed  based  on  a  stress 
analysis  that  accounts  for  the  finite  plate  dimensions  (FIGEOM 
computer  code) .  The  transverse  fastener  displacement  is  governed  by 
a  fourth  order  ordinary  differential  equation,  and  is  computed  by 
solving  a  central  difference  formulation  of  the  equation.  The 
fastener  length  is  divided  into  many  nodes,  each  node  representing  a 
ply  (or  a  fraction  thereof)  in  a  laminate  (see  Figures  10  and  11) . 

The  relative  displacement  between  the  bolted  plates  is 
computed  to  obtain  the  effective  fastener  stiffness  value  as  follows 
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Figure  8.  Boundary  and  Continuity  Conditions  for  a 
Typical  Single-Lap  Joint. 
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Figure  9,  Boundary  and  Continuity  Conditions  for  a 
Typical  Double-Lap  Joint. 
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Figure  10.  An  Example  of  the  Node  Layout  and  Number  Scheme 
in  a  Single  Lap  Shear  Joint  Configuration. 
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Figure  11.  A  General  Node  Arrangement  with  n  Nodes  in  Plate  1 
and  m  Nodes  in  Plate  2. 
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where  uT  and  nB  are  the  average  relative  displacements  between  the 
top  and  bottom  surfaces  in  the  top  and  bottom  plates ,  respectively. 

P  is  the  load  applied  along  the  x  direction  in  the  xy  plane  of  the 
bolted  plates.  When  both  the  bolted  plates  are  isotropic,  the 
affective  transverse  fastener  stiffnesses  in  the  x  and  y  directions 
are  identioal;  i.e.,  ky  ■  kx.  If  either  plate  is  a  laminated 
composite,  this  analysis  is  performed  twice  for  eaoh  fastener.  The 
layup  used  to  compute  kx  is  rotated  by  90  degrees  to  obtain  ky. 

As  shown  in  Figure  12,  tho  two-node  effective  fastener 
element  allows  two  degrees  of  freedom  (DOF)  at  eaoh  node, 
perpendicular  to  the  axis  of  the  element.  Through  the  FDFA 
computation  of  kx  and  ky,  a  4  x  4  effective  fastener  stiffness 
matrix  is  generated.  Until  now,  these  stiffnesses,  referred  to  as 
joint  stiffnesses,  and  were  only  obtained  as  experimentally  measured 
quantities. 

2.2.2  Element  Stiffness  Matrix  For  a  Plate  With  a  Loaded  Hole 

In  a  general  multiply-fastened  panel,  significant  moments 
and  out-of-plane  forces  can  be  generated  by  the  applied  loading, 
particularly  in  a  single- shear  load  transfer  configuration.  The 
present  analysis  assumes  that  the  loaded  hole,  unleaded  hole  and 
plain  elements  behave  essentially  as  membranes  under  plane  stress 
conditions. 

The  characteristic  feature  of  the  loaded  and  unloaded  hole 
elements  is  the  presence  of  a  stross  concentrator  (the  hole)  which 
complicates  the  process  of  determining  stiffness  coefficients.  The 
FIGEOM  code,  developed  and  described  in  Reference  1,  is  capable  of 
computing  the  state  of  stress  within  a  doubly-connected  region  of 
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finite  dimensions,  under  arbitrary  inplane  biaxial  loading.  The 
suooessful  development  of  FIQEOM  motivated  the  adoption  of  a 
flexibility  approach  to  computing  the  element  stiffness  matrix.  The 
natural  mode  method,  originally  proposed  by  Argyris,  is  employed  for 
this  purpose  (Reference  4). 

The  natural  mode  method  was  originally  developed  as  a 
simpler  alternative  to  the  sometimes  tedious  matrix  displacement 
method  of  determining  element  stiffness  relationships.  The  natural 
mode  method  recognises  that  the  total  number  of  kinematic  degrees  of 
freedom  in  an  element  can  be  separated  into  straining  and  rigid  body 
modes.  Only  the  straining  modes  give  rise  to  stiffnesses  that  are 
referred  to  as  natural  or  invariant  stiffnesses.  The  natural 
stiffness  matrix  is  of  a  lover  order  than  the  global  stiffness 
matrix.  The  natural  mode  technique  proceeds  from  a  flexibility 
standpoint  in  which  natural  load  oases  are  initially  imposed  to 
compute  the  natural  flexibilities.  The  natural  flexibility  matrix 
is  subsequently  inverted  to  yield  the  natural  stiffness  matrix.  The 
natural  stiffness  matrix  is  then  expanded  to  yield  the  global 
stiffness  matrix  using  relationships  between  the  natural  modes  and 
the  nodal  displacements. 

The  natural  flexibility  coefficients  are  oomputed  based  on 
the  principle  of  virtual  work,  When  stresses  are  varied  while 
strains  are  held  constant,  a  calculus  of  variations  definition  of 
the  virtual  work  is: 
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where  v  is  the  volume  of  the  domain  of  interest.  The  stresses  and 
strains  introduced  by  the  natural  loads  are  defined  as: 


& 


\>\>v v 


kdl  Jk.  mS  I 


{a}  -  fo  o  T  1T  -  [o  ]  {*»  } 
L  x  y  xyJ  1  J  N 


U>  ■  tcx  ey  ■  hM'1  {V 


whsra  { o }  and  { t  )  ar«  tha  statas  of  atraaa  and  strain  at  a  point  In 
tha  plata  of  thloknass  h,  (pn)  Is  a  vaotor  of  natural  or  ganaralizad 
loads ,  h[A]  Is  tha  inplana  flaxibility  matrix  for  a  laninatad  or 
natallio  plata,  and  <5  )  contains  tha  contribution  of  aach  natural 
load  oasa  to  tha  total  atrass  stata  in  tha  plata.  squation  2  may 
than  ba  vrittan  as: 
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Intagrating  in  tha  thloknass  diraction, 
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than  ba  dafinad  as: 
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where  S  is  the  area  of  the  domain  of  interest*  If  (oN)  is  the 
natural  displacement  vector,  the  flexibility  relationship  is  ex¬ 
pressed  as: 


{V  "  W1  {PN}  “  W  <PN} 


where  [kN  ]  is  the  natural  stiffness  matrix. 


To  relate  the  displacements  in  the  natural  and  global 
coordinate  systems,  the  global  displacement  vector  can  be 
represented  as  a  combination  of  elastic  and  rigid  body  components. 
Assuming  n  nodes  in  the  plate  element,  and  two  degrees  of  freedom  (u 
and  v  in  the  x  and  y  directions,  respectively)  at  each  node, 
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The  elastic  global  displacements  at  the  n  nodes  are  related  to  the 


natural  loads  as  follows: 
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where  [AN  ]  Is  a  transformation  matrix.  Substituting  Equation  9 
into  Equation  11,  ons  obtains: 
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Th*  rigid  body  components  of  th«  global  di.plac.s.nt.  aro 

expressed  as: 


tpc)  -  [*„]  {»;>  (13) 

where  {Pj}  containe  the  rigid  body  translations  in  the  x  and  y 
directions  (u  ,  v  ) ,  and  the  rigid  body  rotation  about  the  z  direc¬ 
tion  (ez).  The  [A0]  matrix  is  solely  dependent  on  the  element 
geometry,  and  an  example  for  a  five-node  element  is  presented  in 
Figure  13. 

The  relationships  in  Equations  12  and  13  are  adjoined  to 
yield  the  following  expression  for  the  global  displacements  (see 
Equation  10) : 


The  invars*  of  Equation  14  yialds  a  ralationship  betwaan  tha 
displacements  in  tha  natural  and  global  coordinata  systems j 
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(15) 
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Incorporating  Equation  16  into  the  principle  of  virtual  work,  tha 
following  relationship  between  the  nodal  (global)  loads  and  the 
natural  loads  is  obtained: 


{P} 


(17) 
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The  global  stiffness  matrix  is  then  related  tc  the  natural  stiff¬ 
ness  matrix  [KN]  through  the  transformation  matrix  [ae],  as  follows: 


(18) 
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The  order  of  the  natural  flexibility  matrix  is  less  than 
the  total  number  of  degrees  of  freedom  (DOF)  in  the  element  by 
three.  The  5-node,  10-DOF  loaded  hole  element,  (See  Figure  14) 
therefore,  requires  seven  natural  load  cases  that  form  an  uncoupled, 
orthogonal  set.  These  load  cases  fully  interrogate  nodal 
interactions,  and  represent  the  basic  element  deformation  modes, 
including  membrane  stretching,  shear  and  banding  (see  Figure  15) . 

In  computing  the  natural  flexibility  matrix,  Equation  7  is  evaluated 
numerically  using  a  standard  Gaussian  integration  scheme  to  approxi¬ 
mate  the  surface  integral: 
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Figure  14.  Five-Node.  Loaded  Hole  Element  With  Depicted  Nodal 
Degrees  of  Freedom 


where  Wj,  and  Wj  are  weighting  functions  at  (x  ,  y  )  locations. 


The  integration  (summation)  in  Equation  19  is  performed  by 
dividing  the  element  into  four  regions.  The  stresses  for  each  load 
case  are  computed  in  each  region,  at  locations  that  correspond  to 
fifth  order  Gaussian  quadrature  points,  scaled  to  the  geometry  of 
the  element.  The  computed  stresses  are  summed  and  weighted  in 
accordance  with  Equation  19  to  yield  the  natural  flexibility 
coefficients.  A  typical  arrangement  of  Gaussian  quadrature  points 
in  an  element  is  shown  in  Figure  16. 


In  the  loaded  hole  element,  the  first  four  load  oases,  in 
which  the  externally  applied  load  is  reacted  at  the  boundary,  cause 
a  significant  non-uniform  distortion  of  the  element  edges.  An  exag¬ 
gerated  displacement  profile  for  one  of  these  load  oases  is  shown  in 
Figure  17.  These  straining  modes  are  not  adequately  represented  by 
storing  only  the  nodal  displacements  in  the  [ANJ  matrix.  To  correct 
this  problem,  the  average  edge  normal  displacements  are  assigned  to 
the  nodes  (see  Figure  17). 


2.2.3  Stiffness  Matrices  For  a  Plate  With  an  Unloaded  (Open) 

Hole  and  For  Plain  (Unnotched)  Elements 

The  generation  of  global  stiffness  matrices  for  the  open 
hole  and  plain  elements  follows  the  procedure  outlined  in  Section 
2.2.2.  These  elements  contain  only  four  nodes  (8  DOF)  each  (see 
Figure  18) .  Therefore,  only  five  natural  load  cases  are  required  to 
generate  their  natural  flexibility  matrices  (Figure  19)  .  The  trans¬ 
formation  of  the  5x5  natural  flexibility  matrices  to  the  8x8 
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Figure  18. 


The  Four-Node  Open  Hole  Element  With  Depicted  Nodal  Degrees  of 
Freedom 
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global  stiffness  matrices  for  the  two  elements  follows  Equations  9 
to  18.  Sinoa  FIGEOM  was  developed  to  analyze  doubly-connected 
planform  regions,  the  plain  element  stiffnesses  are  obtained  using 
the  open  hole  element  algorithm,  setting  the  hole  radius  to  a  mini¬ 
mum  value.  This  also  provides  the  added  benefit  of  using  the  same 
set  of  subroutines  to  generate  the  stiffness  matrices  for  all  the 
plate  elements  (loaded  hole,  unloaded  hole  and  plain  elements) . 

2.3  Load  Distribution  Among  Fasteners 

Figure  20  shows  a  typical  bolted  joint  and  its  finite 
element  representation.  The  stiffnesses  of  each  element  are 
initially  computed  and  stored  (see  Section  2.2).  The  global  joint 
stiffness  matrix  is  then  formed  by  assembling  the  individual 
stiffness  matrices  for  the  loaded  hole,  unloaded  hole,  plain  and 
effective  fastener  elements.  As  described  in  Section  3.1,  the  SAMCJ 
user  only  defines  the  type,  geometry  and  properties  of  the 
individual  elements.  SAMCJ  internally  processes  this  information  to 
genorate  the  global  joint  stiffness  matrix. 

Figure  3  shows  the  load  introduction  and  boundary  con¬ 
straint  locations  assumed  by  SAMCJ.  Currently,  SAMCJ  assumes  that  a 
uniaxial  l-kip  load  is  applied  along  the  x  (longitudinal)  direction 
on  the  left  edge  of  the  top  plate  (see  Figure  3).  Therefore,  the  u 
displacement  is  constrained  along  the  right  edge  of  the  bottom 
plate.  A  cornar  node  along  this  edge  is  also  constrained  in  the  y 
direction  to  prevent  a  rigid  body  translation  in  that  direction  (see 
Figure  3).  The  user  specifies  the  applied  load  to  be  tensile  or 
compressive. 

The  assembled,  global  joint  stiffness  matrix  is  related  to 
the  nodal  displacements  and  loads  as  follows: 


(!>}  »  [K  1  (p) 

8  (20) 
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Element  Model  of  a  Sample  Taoered  Bolted  Joint 


Each  bolted  plate  contains  M  nodes,  which  include  the  fastener 
nodes.  Each  node  has  two  degrees  of  freedom  (u  and  v 
displacements) .  Therefore,  the  global  joint  stiffness  matrix  is  4M 
x  4M  in  sizo.  The  imposed  boundary  constraints  reduce  the  size  of 
the  stiffness  matrix  that  is  eventually  used  to  compute  the  nodal 
displacements  and  boundary  constraint  forces.  Incorporation  of  the 
nodal  displacements  into  element  equilibrium  equations  yield  the 
nodal  forces  in  the  individual  elements.  The  loads  at  node  5  in 
eaah  loaded  hole  element  provide  the  x  and  y  components  of  the 
fastener  load  in  that  element. 

2.4  Stress  State  at  Any  Location  in  a  Bolted  Plate 

The  stress  stats  at  any  internal  point  in  an  element  is 
computed  using  a  procedure  similar  to  the  computation  of  the  natural 
flexibility  matrix.  During  the  generation  of  element  natural  flsxi- 
bility  matrices,  the  stress  states  at  locations  within  the  element  are 
computed  and  stored  for  unit  values  of  every  natural  load  case.  Desired 
stress  recovery  locations  are  pre-selected  for  this  purpose.  The 
relationship  between  the  stress  states  at  these  pre-selected  n  points 
and  the  natural  force  system  is: 


(21) 


where  contains  a x  ,  ov  and  oxy  at  the  selected  S  locations,  and 
[sN]  contains  the  stresses  per  unit  natural  load  at  the  same 
locations.  As  discussed  in  Section  2.2.2,  the  natural  load  vector  is 
related  to  the  natural  nodal  displacement  vector  as  follows: 
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And,  the  natural  displacements  are  related  to  the  global  (nodal) 
displacement  as  follows: 


(PN>  ■  [*e3T{p} 


(23) 


where  (24) 

{p}  -  (ux  Vl  u2  v2  .  .  .  .  u10  v1Q} 

Substituting  Equations  21  and  22  into  Equation  20,  one  obtains: 


to} 


[sj  cs,]  0.]T 
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(25) 


where  [s]  ie  a  3N  x  10  matrix  for  a  loaded!  hole  element  (io  DOF) ,  and 
a  3N  x  8  matrix  for  an  open  hole  element  (8  DOF). 


Terms  in  the  [g]  matrix  are  computed  at  the  same  time  the 
element  flexibility  and  stiffness  matrices  aro  calculated.  Hence, 
after  the  global  joint  equations  are  solved  for  the  nodal  displace-’ 
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merits,  the  stresses  in  any  element  are  recovered  at  these  prese¬ 
lected  points  using  Equation  24. 

2.5  Strength  and  Failure  Mode  Prediction  for  a  Bolted  Plate 

SAMCJ  predicts  the  strength  of  a  bolted  plate  using 
average  stress  failure  criteria  at  the  laminate  level  (see  Figure 
21) .  In  Reference  1,  the  same  criteria  were  applied  at  the  lamina 
level  to  predict  progressive  ply  failures  in  a  singly-fastened 
laminate  (the  SASCJ  computer  code) .  The  characteristic  distances 
over  which  the  stresses  are  averaged  in  SAMCJ  are  different  from 
those  used  in  SASCJ.  Also,  failure  is  assumed  to  be  a  one-step 
(catastrophic)  process.  The  strength  of  a  bolted  plate  corresponds 
to  the  initial  failure  at  a  fastener  or  cut-out  location,  in  the 
bearing,  shear-out  or  net  section  failure  mode. 

Appropriate  failure  sites  are  identified  by  SAMCJ  in  every 
element.  These  sites  vary  with  the  applied  loading  (see  Figure  22). 
The  characteristic  distances  for  the  net  section,  shear-out  and 
bearing  modes  of  failure  are  divided  into  many  regions.  Following 
the  procedure  outlined  in  Section  2.4,  the  appropriate  stress 
components  corresponding  to  a  1-kip  joint  load  are  computed  at  these 
points,  and  their  average  values  over  the  respective  characteristic 
distances  are  stored.  The  ratios  of  these  average  stresses  to  the 
corresponding  unnotched  strengths  are  subsequently  computed  and 
relatively  evaluated  to  predict  the  strength  of  the  bolted  plate, 
the  failure  site  and  the  failure  mode.  Under  tensile  loading,  the 
average  ~x  value  over  ans  is  divided  by  the  unnotched  tensile 
strength  to  predict  net  section  tensile  failure.  Under  compressive 
loading,  the  avarage  cx  value  over  ans  is  divided  by  the  unnotched 
compressive  strength  to  predict  net  section  compressive  failure. 

The  average  ;x  value  over  abr®  is  divided  by  the  unnotched  compres¬ 
sive  strengch  to  predict  bearing  failure.  The  average  ;xy  value 
over  aso  is  divided  by  the  unnotched  shear  strength  to  predict 
shear-out  failure. 
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Figure  21  .  The  Characteristic  Distances  Used  in  the  Average  Stress 


Failure  Criteria. 


The  unnotched  laminate  strengths,  under  tension  and  under 
compression,  are  computed  by  SAMCJ  based  on  input  fiber-directional 
failure  strain  values  (tensile  and  compressive) .  Laminate  strengths 
under  Nx  and  Nxy  loadings  (inplane  normal  and  shear  stress 
resultants,  respectively)  are  assumed  to  correspond  to  first  fiber 
failure  in  a  ply.  This  simplistic  strength  prediction  procedure 
introduces  inaccuracies  that  have  been  acknowledged  and  discussed  in 
the  literature.  Nevertheless,  SAMCJ  adopts  this  procedure  for  lack 
of  a  validated  alternative. 


SAMCJ  assumes  that  a  net  section,  shear-out  or  bearing 
failure  of  any  element  results  in  joint  failure.  This  assumption 
results  in  a  one-step  strength,  failure  site  and  failure  mode  pre¬ 
diction  for  a  multiply-fastened  plate.  A  1-kip  tensile  or  compres¬ 
sive  load  is  applied,  and  fastener  loads  and  normalized  averaged 
stresses  corresponding  to  net  section,  shear-out  and  bearing  failure 
modes  are  computed.  The  failure  value  of  the  applied  load  corres¬ 
ponds  to  a  unit  value  of  the  maximum  normalized  average  stress.  An 
identification  of  the  maximum  normalized  average  stress,  and  its 
location,  provides  the  joint  failure  mode  and  the  critical  fastener 
or  cut-out  location. 

2.6  current  SAMCJ  Limitations 

SAMCJ  is  a  versatile  code  that  predicts  the  strength  of 
bolted  laminates  and  the  corresponding  failure  mode.  However,  there 
are  segments  of  the  analysis  that  can  bo  improved  through  additional 
efforts  that  are  beyond  the  scope  of  this  program.  These 
limitations  are  discussed  below. 

The  representation  of  a  loaded  hole  by  a  five-node  element 
(with  ten  DOF)  is  discussed  in  Section  2.2.2.  Let  a  and  b  be  the 
planform  dimensions  of  the  element  in  the  x  and  y  directions, 
respectively,  and  D,  the  hole  diameter.  When  a/D  and  b/D  are  small 
(less  than  2),  and  the  element  aspect  ratio  (a/b)  is  less  than 


unity,  the  highly  distorted  element  shape  cannot  be  appropriately 
accounted  for  even  with  averaged  displacements  assigned  to  the  nodes 
(see  Figure  17) . 

Fastener  load  distributions  were  predicted  for  two  test 
cases  using  three  forms  of  the  [A0]  matrix.  In  one  form,  the  actual 
nodal  displacements  were  used  as  the  A  terms.  In  the  second  form, 
the  midside  displacement  value  for  the  highly  deformed  edge  was 
assumed  to  be  the  nodal  values  on  that  side  in  the  [A0]  matrix.  In 
the  third  form,  the  average  displacement  of  the  deformed  edge  was 
used  ns  the  nodal  values  in  the  [A0]  matrix.  The  considered  test 
cases  are  from  Reference  2.  One  contains  two  fasteners  in  the  load 
direction,  and  the  other  contains  five  fasteners  in  the  load 
direction.  Figures  23  and  24  present  the  fastener  load 
distributions  predicted  by  the  three  forms  of  the  [A0]  matrix,  along 
with  experimentally  measured  fastener  loads  for  the  two  test  oases 
(Reference  2) .  The  best  correlation  between  analytical  predictions 
and  test  measurements  is  obtained  when  the  nodal  displacements  are 
assumed  to  be  the  average  value  along  the  deformed  edge.  SAMCJ, 
therefore,  generates  the  [A0]  matrix  for  the  loaded  hole  element 
following  this  procedure. 

In  dividing  a  bolted  plate  into  many  elements  (loaded  or 
unloaded  hole  elements,  as  well  as  plain  elements) ,  it  is  advisible 
to  maintain  element  geometries  that  do  not  render  the  generated 
stiffness  matrices  inaccurate.  Figure  25  presents  results  from  a 
study  conducted  on  a  singly  fastened  metallic  plate.  P  is  the 
recovered  load  that  is  obtained  by  integrating  the  stress  along  a 
line  transverse  to  the  load  direction  as  shwon  in  Figure  25.  P 
is  the  applied  load,  or  the  sum  of  the  nodal  loads  (especially  in 
the  interior  elements  in  a  general  multifastened  plate) .  The 
recovered  load  (P  )  approaches  the  applied  load  value  (P)  when  the 
plate  aspect  ratio  (a/b)  increases  beyond  unity.  Also,  a/D  and  b/D 
must  have  a  minimum  value  of  approximately  3.  In  predicting 
failure  in  the  net  section,  bearing  and  shear-out  modes,  the 
computed  average  stress  values  are  multiplied  by  P/Pr,  to  remove 
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geometry  (modeling)  effects  from  the  computod  stresses. 

The  generetion  of  e  higher  order  loaded  hole  element  may 
eliminate  the  approximation  introduced  by  the  five  node  element,  and 
result  in  a  Pr/P  value  that  ie  approximately  unity  for  any  element 
geometry.  A  nine  node  element,  including  midside  nodes,  is 
recommended  for  future  investigation.  Other  factors  that  will 
improve  SAMCJ  predictions  are  more  accurate  computations  of  the 
unnotched  strengths,  and  a  modified  fastener  analysis  that  can 
account  for  oountersunk  fastener  geometry. 

2.7  Design  Application 

The  design  of  a  bolted  joint  in  composite  structures 
involves  the  selection  of  the  fastener  type,  size  and  arrangement 
(spaaing  between  adjacent  fasteners),  and  geometry  changes  in  the 
bolted  plates  (layup  change  and  change  in  the  planform  dimensions). 
SAMCJ  can  quickly  interrogate  the  effects  of  all  these  parameters  on 
the  joint  strength,  to  provide  a  near  optimum  value  for  eaoh.  In 
doing  so,  SAMCJ  is  independent  of  test  measurements  like  "joint 
stiffnesses,"  and  is,  thsrefore,  a  rapidly  usable  analytical  design 
tool.  If  the  bolted  laminate  is  to  be  fabricated  using  a  new 
material,  only  the  basic  lamina  properties  and  the  characteristic 
distances  for  the  average  stress  failure  criteria  have  to  be 
determined,  prior  to  performing  the  analysis.  When  a  charaoterized 
material  is  used  in  the  boltsd  structural  part,  SAMCJ  predicts  its 
strength  without  requiring  complementary  test  rosults.  The  fastener 
size  and  spaaing,  and  the  boltsd  plate  geometry  and  properties,  are 
varied  systematically  to  analytically  predict  their  effect  on  the 
joint  strength  and  efficiency  (weight,  durability,  etc.).  These 
parametric  studies  provide  guidance  in  the  selection  of  the  most 
efficient  bolted  joint  configuration  for  the  assumed  loading. 

2 . 8  Test  Requirements 

As  mentioned  in  Section  2.7,  SAMCJ  predicts  the  strengths 
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of  bolted  laminates  without  requiring  complementary  test  results 
when  the  laminate  is  fabricated  using  a  characterised  material.  A 
material  is  said  to  be  characterised  when  the  basic  lamina 
properties  (stiffnesses,  strengths,  failure  strains  and  other 
physiochemomeohanical  properties)  and  relevant  structural  properties 
(characteristic  distances  for  the  prediction  of  net  section,  bearing 
and  shear-out  modes  of  failure  in  notched  laminates,  eto.)  are 
available.  Zf  a  new  material  is  used  in  a  bolted  laminate,  its 
basic  lan in a  properties  and  failure  parameters  have  to  be  obtained 
prior  to  using  8AMCJ  for  strength  predictions. 
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SECTION  3 


ANALYTICAL  PREDICTIONS 

In  ths  experimental  part  of  Tack  2  in  this  Northrop/AFWAL 
program,  ovar  160  oompoa ite-to-aluminum  multifaatanar  joints  wars 
tested  undar  atatio  loading  (aaa  Tabla  l  and  Figure  26) .  Rafaranoa 
2  contains  raaults  from  thaas  tests,  including  faatanar  load 
distribution  maaauramanta  uaing  strain-gaged  bolts,  failura  loads, 
failura  locations  and  failuro  modes.  Sampla  test  oasas  from 
Referanca  2  ara  analyzed  below  using  tha  SAMCJ  computer  ooda. 
Analytical  predictions  ara  compared  with  test  results  to  establish 
tha  validity  of  the  developed  analysis. 

Bolted  laminates  ware  fabricated  using  ASl/3501-6 
graphita/apoxy  unidirectional  praprag  material  containing 
approximately  35%  resin  by  weight.  These  included  20-  and  40-ply 
laminates  with  50/40/10,  70/20/10,  30/60/10  and  25/60/15 
(percentages  of  0  ,  +-45  and  90  plies,  respectively)  layups. 

The  20-p.ly  30/40/10,  70/20/10  and  30/60/10  layups  had 
[45/0/-45/0) 2/0/90]  ,  [ 45/0/-45/03/90/03 ] ■  and  [ 45/0/-45/0/45/90/ 

-4 5/0/+-45 ] s  stacking  sequence*,  respectively.  The  40-ply  50/40/10 
and  70/20/10  layups  had  [45/0/-45/0) 2/0/90] 2s  and 

[4 5/0/-45/03/90/03 ] 2s  stacking  sequences,  respectively.  The  40-ply 
25/60/15  laminate  had  a  [45/0/-45/0/45/90/-45/0/+-45] 2s  stacking 
sequence,  with  the  twelfth  0  ply  replaced  by  a  90  ply. 

The  tested  fantener  arrangements  included:  two  fasteners 
in  tandem,  two  at  an  angle  to  the  load  direction,  three  fasteners  in 
tandem,  three  fasteners  in  each  of  two  columns  with  an  adjacent  cut¬ 
out,  and  four  fasteners  in  each  of  two  columns  with  a  cut-out  either 
between  or  adjacent  to  the  rows.  Rows  and  columns  of  faoteners  are 
along  and  perpendicular  to  the  load  direction,  respectively.  The 
fastener  spacing  in  the  load  and  transverse  directions  (SL  and  ST, 
respectively) ,  specimen  width  and  edge  distance  (W  and  E, 


11  ti'sr  cases  Duxeter  5/16  inch  (PU);  Torque  =  100  in-lb;  Test  Eaylronaait  is  KTP;  and  Joint  Type  I*  Single-Lap  (SL)  unless  otherwise  noted 

Static  Tension  PH  «  Protruding  Head  SL  ■  Single— Lap  KTD  *  Icoa  T«peratu re  Dry 

Static  Compression  CSV.  -  100°  Countersink  DL  *  Double-Lap  ETW  *  Elevated  T«perature  Wet 

Tension  Head 


respectively) ,  and  cut  out  diameter  (H  )  and  location,  for  the 
various  test  cases,  are  listed  in  Table  1. 

A1  minum  plates  were  bolted  to  laminates  to  effect  load 
transfer  in  single-  and  double-shear  configurations.  The  metallic 

plates  were  machined  from  7075-T7  raw  stock,  and  contained  fastener 
hole  arrangements  that  were  compatible  with  those  in  the  laminated 
specimens.  Figure  26  presents  the  dimensions  of  the  metal  plates 
used  in  the  various  tests. 

Most  of  the  tests  used  5/16-inch  diameter,  protruding 
head  steel  fasteners.  Selected  tests  in  a  single  shear 
configuration  used  5/16-inch  diameter,  100  countersunk  (tension 
head)  steel  fasteners.  The  fasteners  were  torqued  to  100  in. -lb, 
prior  to  testing,  unless  otherwise  specified. 

3.1  Composite-to-Metal  Joints  with  Two  Fasteners  in  Tandem 

Four  test  cases  (201,  202,  203  and  206  in  Table  1), 
addressing  joints  with  two  fasteners  in  tandem  (along  the  loading 
direction),  were  analyzed  using  the  SAMCJ  computer  code.  Test  case 
201  was  modeled  as  shown  in  Figure  27.  The  top  plate  was  modeled  to 
be  the  0.31-inch-thick  aluminum  plate,  and  the  bottom  plate  was 
modeled  to  be  the  20-ply,  50/40/10  (percentages  of  0  /+-45  /90 
plies)  graphite/epox.y  plate.  A  1 -kip  load  is  applied  to  the  left 
edge  of  the  top  plate,  and  the  right  edge  of  the  bottom  plate  is 
constrained,  as  shown  in  Figure  27.  u  and  v  are  the  nodal 
displacements  in  the  x  and  y  directions,  respectively.  The 
mentioned  loading  of  the  top  plate  and  the  constraining  of  the 
bottom  plate  are  automatically  done  by  SAMCJ.  If  two  or  more  rows 
of  elements  were  present  in  the  model,  the  applied  1-kip  load  will 
be  distributed  among  three  or  more  nodes  in  proportion  to  the 
element  widths. 

The  SAMCJ  input  data,  in  English  units,  for  test  case  201 
are  listed  in  Figure  28.  The  first  entry  (1)  identifies  the  loading 
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*  GRID-POINT  IDENTIFICATION 
*>'<  LOADED  HOLE  ELEMENT  IDENTIFICATION 
+  EFFECTIVE  FASTENER  ELEMENT  IDENTIFICATION 


BOTTOM  PLATE  (GRAPHITE/EPOXY) 
204 


1.875 


0.9375 


** 

401 

y 

402 

2(17__0 

\jLx  208— O 

501 

0°  DIRECTION  502 

201  i 

jsn 

.  205 

0.9375  1.25 


2.1873  2.5 


Floury  27.  Finite  Element  Model  for  Test  Case  201 
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Figure  28.  SAMCJ  Input  for  Test  Case  201 


configuration  to  be  a  single-shear  configuration.  For  a  double¬ 
shear  load  transfer,  this  entry  would  be  2.  The  second  entry  (1) 
identif  <«  the  load  to  be  in  static  tension.  For  static  compressive 
loading,  unis  entry  would  be  2.  The  third  entry  requests  the  type 
of  fastener  used  in  the  joint.  The  entry  of  a  1  specifies  the 
fastener  to  have  a  protruding  head  and  a  2  specifies  the  fastener  to 
have  a  countersunk  head.  The  next  two  entries  say  that  the  top 
plate  is  a  metal  (M) ,  identified  as  "Aluminum."  The  two  entries 
following  these  say  that  the  bottom  plate  is  a  composite  laminate 
(C) ,  identified  as  "50/40/10  ASl/3501-6. . . ".  Subsequently,  the 
Young's  modulus  (10.0D6)  and  Poisson's  ratio  (0.3)  for  aluminum,  and 
the  fiber-direotional,  transverse  and  shear  moduli  and  Poisson's 
ratio  (18.5D6,  1.906,  0.85D6  and  0.3,  respectively)  for  the 
composite  lamina  are  input.  The  next  five  entries  specify  that  four 
(4)  different  fiber  orientations  are  present  in  the  laminate  (0,  45, 
-45  and  90  degrees  with  respect  to  the  loading  direction) .  The 
following  three  entries  say  that  the  elemente  in  the  bottom  plate 
contain  1  layup  of  20  plies,  of  0.565-inch  thickness  each.  The 
stacking  sequence  for  tnis  layup  is  input  next,  where  1,  2,  3 
and  4  refer  to  0,  45,  -45  and  90-degree  fiber  orientations, 
respectively.  Subsequently,  the  fastener  ia  idantifisd  as  "Steal, " 
and  its  Young's  modulus,  Poisson's  ratio  and  diameter  (30.0D6,  0.3, 
0.3125)  are  input. 

Eight  grid-points  each  are  specified  in  the  top  and 

bottom  plates  (101  to  108  and  201  to  208,  respectively),  along  with 
their  x  and  y  coordinates  (see  Figure  27).  Following  this,  two 
elements  are  specified  in  each  plate,  along  with  their  nodal 
connectivity  and  element  type  information.  Nodal  connectivity  is 
specified  starting  from  the  bottom  left  node,  going  clockwise  around 
the  element  boundary,  and  ending  at  the  fastener  (internal)  node. 
Element  301  in  the  top  plate,  for  example,  has  101,  102,  104  and  103 
as  ite  corner  nodes,  and  107  as  its  fastener  node.  The  fifth  node 
will  be  entered  as  0  for  plain  and  unloaded  hole  elements.  The 
element  type  information  follows  the  fifth  node  identification.  It 
ia  1,  2  and  3  for  plain,  loaded  hole  and  unloaded  hole  elements, 


respectively.  Following  this,  additional  element  data  are  specified 
for  the  two  plates.  These  include  the  element  thicknesses  or  layup 
identification  number  for  plain  and  loaded  hole  elements,  and 
additional  information  (x  and  y  coordinates  of  the  hole  center 
and  the  hole  radius)  for  unloaded  hole  elements.  For  test  case  201, 
elements  301  and  302  in  the  top  plate  (metal)  are  specified  to  be 
0.3l-inch  thick.  Elements  401  and  402  in  the  bottom  plate  (composite) 
are  specified  to  contain  the  stacking  sequence  identified  as  one  (1). 
The  element  definitions  are  succeeded  by  the  definition  of  two 
effective  fasteners  (501  and  502) .  These  are  identified  as  fasteners 
that  connect  nodes  107  and  108  in  the  top  plate  to  nodes  207  and 
208  in  the  bottom  plate,  respectively. 

The  one  (1)  following  this  states  that  groups  of  identical 
elements  will  be  specified  in  the  two  plates.  If  two  (2)  is  entered 
here,  all  elements  will  be  assumed  to  be  different  from  one  another, 
resulting  in  larger  computational  costs.  The  entry  "110  0"  refers 
to  the  number  of  groups  of  effective  fasteners,  loaded  hole,  unloaded 
hole  and  plain  elements,  respectively,  in  the  top  plate.  A  zero  (0) 
specifies  the  absence  of  an  element  type.  The  number  of  elements  in 
each  group,  and  the  corresponding  element  numbers,  are  input 
subsequently.  In  Figure  27,  two  identical  loaded  hole  elements 
(301  and  302)  are  identified  in  the  top  plate.  Following  this,  the 
number  of  groups  of  loaded  hole,  unloaded  hole  and  plain  elements 
in  the  bottom  plate  (1,  0  and  0,  respectively)  is  entered.  In  figure 
27,  2  identical  loaded  hole  elements  (401  and  402)  are  identified 
in  the  bottom  plate. 

The  last  four  lines  of  input  introduce  the  failure 
parameters  for  the  materials  in  the  two  plates.  For  metallic 
plates,  the  tensile,  compressive  and  shear  strengths  (250. 0D3  each 
in  Figure  27) ,  and  the  averaging  distances  for  net  section,  bearing 
and  shear-out  modes  of  failure  (0.5  each  in  Figure  28)  are  input, 
since  the  joints  were  designed  to  fail  the  laminated  plates,  and 
SAMCJ  was  developed  primarily  for  the  prediction  of  the  strength  of 
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bolted  laminates,  the  failure  parameters  for  the  metallic  plates 
were  input  to  be  arbitrarily  high.  This  information  is  followed  by 
the  failure  parameters  for  the  bottom  (composite)  plate.  The  first 
line  specifies  the  fiber  directional  failure  strains  for  the  material 
under  tension  (0.012),  compression  (0.0175)  and  shear  (0.012).  These 
values  are  used  by  SAMCJ  to  compute  the  plain  laminate  tensile, 
compressive  and  shear  strengths,  based  on  laminated  plate  theory  and 
the  assumption  of  laminate  failure  corresponding  to  first  fiber 
failure  in  any  of  ita  plies.  Ths  last  lins  in  Figure  28  epeoifiee 
the  distances  over  which  the  longitudinal  (0.10  and  0.25)  and  shear 
(0.25)  stress  components  ars  avaraged,  to  prediot  net  section, 
bearing  and  shear-out  modes  of  failurs,  respectively  (see  Figure 
21)  . 


SAMCJ  predictions  for  test  case  201  are  compared  with  test 
results  from  Reference  2  in  Figure  29.  The  two  fasteners  are 
predicted  to  carry  nearly  equal  loads,  in  agreement  with  the  values 
measured  using  atrain-gaged  bolte  (Reference  2) .  The  assumed 
failure  parameters  (aons ,  a0brg  and  aoso  of  0.10,  0.25  and  0.25 
inch,  respectively)  and  the  average  stress  failure  criteria  predict 
the  observed  shear-out  mode  of  failure  at  the  inner  fastener 
location(l).  The  predicted  failure  load  (9.48  kips)  agrees  very 
well  with  the  average  measured  value  (9.51  kips).  Note  that  a 
bearing  mode  of  failure  at  the  same  fastener  location  is  predicted 
to  occur  at  only  a  slightly  larger  load  level  (9.71  kips).  This 
indicates  the  possibility  of  either  failure  mode,  within  the  scatter 
region  of  the  measured  failure  load.  Observations  in  Reference  2 
also  indicated  the  dual  failure  mode  possibility  for  specimens  in 
test  case  201. 

SAMCJ  predictions  for  test  case  202  are  compared  with  test 
results  in  Figure  30.  In  this  case,  the  model  in  Figure  27  is  used 
along  with  a  70/20/10  layup  for  the  bottom  plate.  The  input  data  in 
Figure  28  can  be  easily  modified  to  account  for  this  change.  SAMCJ 
predicts  a  nearly  equal  load  distribution  between  the  two  fasteners, 
and  a  shear-out  mode  of  failure  in  the  highly  fiber-dominated  layup. 
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Figure  2'i.  SAMCJ  Predict ions  and  Test  Results  for  Test  Case  201 
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Figure  30.  SAMC.I  Predictions  and  Test  Results  for  Test  Case  202 
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As  in  tost  case  201,  and  in  every  other  test  case  discussed  in  this 
section,  the  failure  parameters  were  assumed  to  be  invariant  (a 
values  of  0.10,  0.25  and  0.25  inch  for  net  seotion,  bearing  and 
shear-out  modes  of  failure,  respectively).  The  predicted  failure 
load  (6.90  kips)  is  131  lower  than  the  average  measured  value  (7.92 
kips) . 


Figure  31  presents  SAMCJ  predictions  for  test  case  203, 
where  the  laminate  layup  is  changed  to  a  30/60/10  layup.  The  input 
data  for  this  case  are  identical  to  those  for  the  previous  two 
cases,  with  the  exception  of  the  laminate  layup  definition  and  its 
cured  ply  thickness.  In  this  case,  a  bearing  mode  of  failure  is 
predicted  at  the  inner  fastener  location.  This  agrees  with  the 
observation  in  Reference  2.  The  predicted  failure  load  for  this 
case  (7.18  kips)  is  20%  lower  than  the  average  measured  value  (8.92 
kips) . 


The  compos ite-to-metal  joint  in  test,  case  206  is  tested 
in  a  double  shear  configuration,  instead  of  the  single-shear 
configuration  in  test  case  201.  Test  results  from  Reference  2  are 
compared  with  SAMCJ  predictions  for  this  case  in  Figure  32.  Again, 
the  measured  nearly  equal  load  distribution  is  predicted  by  SAMCJ, 
along  with  the  observed  shearout  mode  of  failure  at  the  inner 
fastener  location.  The  predicted  failure  load  (9.57  kips)  is  in 
excellent  agreement  with  the  measured  average  value  (9.60  kips). 

3.2  Compos ite-to-Metal  Joints  with  Two  Fasteners  at  an  Angle 

to  the  Load  Direction 

The  bolted  plates  in  test  cases  216  to  224  (see  Table  1) 
contain  two  fasteners  each  at  an  angle  to  the  load  direction,  when 
the  fastener  spacings  in  the  load  direction  (S  )  and  the  transverse, 
direction  (ST)  are  reduced,  the  modeling  of  each  bolted  plate 
results  in  inaccurate  SAMCJ  computations.  Refer  to  Figure  33  for 
three  different  models  of  each  bolted  plate  in  test  case  22.1  (S  /n 
S  /D  -  4) .  The  off-center  fastener  location  in  the  four  element 
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models  of  the  bolted  plate  causes  two  of  the  element  boundaries  to 
be  very  close  to  the  fastener  hole  boundary.  These  distances  are 
reduced  even  further  in  other  test  cases  (216  to  220,  and  222  to  224 
in  Table  1) .  Consequently,  the  computed  average  stresses  at  the 
fastener  locations  are  influenced  by  the  proximity  of  the  element 
boundaries . 


Two-,  four-  and  sixteen-element  models  of  each  bolted  plate 
in  tost  case  221  result  in  SAMCJ  predictions  that  are  shown  in 
Figure  34.  In  agreement  with  the  strain-gaged  bolt  measurements  in 
Reference  2,  the  two  fasteners  are  predicted  to  carry  equal  loads 
that  are  aligned  along  the  load  direction.  While  the  equal  load 
distribution  is  to  be  expected,  the  lack  of  a  significant  transverse 
(y)  component  of  the  fastener  load  defeats  intuition.  SAMCJ 
predictions  and  the  test  results  from  Reference  2  indicate  nearly 
zero  values  for  and  ,  pointing  to  the  absence  of  significant  y 
components  of  the  fastener  loads. 

SAMCJ  predicts  the  failure  location  (inner  fastener)  and 
the  failure  mode  (shearout)  recorded  during  the  test  In  Reference  2. 
However,  the  predicted  failure  load  is  nonconservative,  and  is 
dependent  on  the  modeling  of  the  bolted  plates.  For  the  two-element 
models  in  Figure  33,  SAMCJ  predicts  a  failure  load  of  12.5  kips, 
which  is  29%  larger  than  the  measured  average  value  (9.7  kips).  For 
the  four-  and  sixteen  element  models,  the  SAMCJ  predictions  (15.2  and 
16.6  kips)  are  57%  and  71%  larger  than  the  measured  average  value, 
respectively. 

When  the  fasteners  are  very  close  to  one  another,  as  in 
test  case  216  to  224,  SAMCJ  is  unable  to  conservatively  predict  the 
measured  failure  loads.  A  higher  order  element,  proposed  for  future 
development  in  Section  2,  is  expected  to  yield  better  analytical 
predictions  for  such  test  cases. 

3.3  Composite-to-Metal  Joints  with  Four  Fasteners  in  a 

Rectangular  Pattern 
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Test  Case  221,  Static  Tension,  Single-Shear 
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Test  cases  225,  229  and  230  in  Tablo  1  consider  composite- 
to-metal  joints  in  a  single-shear  configuration,  with  four  fasteners 
in  a  rectangular  pattern.  The  finite  element  model  for  each  of  the 
bolted  plates  in  test  cases  225,  229  and  230  is  shown  in  Figure  35. 
The  laminates  in  the  three  test  cases  contained  20  plies  each  in 
50/40/10,  70/20/10  and  30/60/10  layups. 

Figure  36  compares  SAMCJ  predictions  with  the  test  results 
from  Reference  2  for  test  case  225.  The  predicted  fastener  load 
distribution  is  nearly  equal,  with  the  inner  fasteners  carrying  a 
slightly  larger  fraction  of  the  load  (8%) .  This  agrees  with  the 
test  results  in  Reference  2.  SAMCJ  predicts  a  shearout  failure  at 
the  inner  fastener  location.  The  predicted  shearout  load  levels 
corresponding  to  the  two  inner  fasteners  were  18.1  and  18.9  kips. 

In  Reference  2,  the  same  failure  location  (inner  fasteners)  and 
failure  mode  were  observed  in  two  out  of  three  replicates.  In  the 
remaining  replicate,  though,  a  net  section  failure  occurred  across 
the  inner  fasteners.  The  shearout  failures  in  two  specimens  were 
accompanied  by  delaminations  in  the  laminate.  The  failure  load 
predicted  by  SAMCJ  (18.1  kips)  is  only  6%  larger  than  the  measured 
average  value  (17.1  kips). 

SAMCJ  predictions  for  the  bolted  70/20/10  laminate  in  test 
case  229  are  presented  in  Figure  37.  The  fastener  load  distribution 
is  identical  to  that  predicted  for  the  50/40/10  laminate,  and  agrees 
with  the  test  results  from  Reference  2.  The  predicted  railure 
location  (inner  fasteners)  and  failure  mode  (shearout)  correlate 
well  with  the  observations  in  Reference  2.  As  before,  shearout  is 
accompanied  by  dalaminations  in  the  failed  laminata.  The  failure 
load  predicted  by  SAMCJ  (11.8  kips)  is  21%  lower  than  the  measured 
average  value  (14.9  kips), 

SAMCJ  predictions  for  the  bolted  30/60/10  laminate  in  test 
case  230  are  presented  in  Figure  38.  In  this  case,  the  predicted 
loads  in  the  inner  fasteners  (3,4)  are  17%  larger  than  those  in  the 
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Test  Case  230,  Static  Tension,  Single-Shear 
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outer  fasteners  (1,2).  SAMCJ  predicts  a  net  section  failure  across 
the  inner  fastener  holes  (3,4)  at  a  load  level  of  12.4  kips.  In 
Reference  2,  the  predicted  net  section  failure  was  observed  at  the 
predicted  site,  accompanied  by  delaminations.  The  failure  load 
predicted  by  SAMCJ  (12.4  kips)  is  24%  lower  than  the  measured 
average  value  (16.4  kips). 

3.4  Composite-to-Metal  Joints  with  Three  Fasteners  in  a 

Triangular  Pattern 

Test  cases  233  to  241  in  Table  1  address  double- shear  load 
transfer  joints  that  contain  three  fasteners  in  a  triangular 
pattern.  In  predicting  the  effect  of  these  fastener  patterns,  SAMCJ 
encountered  the  same  difficulties  described  in  Section  3.2.  Figure 
39  presents  a  four-element  model  of  the  bolted  plates  in  test  cases 
234,  238  and  239.  SAMCJ  used  this  model  to  make  the  predictions 
discussed  below. 

A  50/40/10  laminate  is  considered  in  test  case  234.  The 
fastener  loads  and  their  orientations  to  the  load  direction, 
predicted  by  SAMCJ,  are  presented  in  Figure  40.  SAMCJ  predictions 
correlate  well  with  the  strain-gaged  bolt  measurements  in  Reference 
2.  Test  observations  indicated  a  combination  of  shearout  and  net 
section  failures  at  the  innermost  fastener  location.  SAMCJ  predicts 
a  shearout  failure  at  the  same  location.  The  failure  load  predicted 
by  SAMCJ  (17.0  kips)  is  38%  larger  than  the  measured  average  values 
(12.3  kips).  As  in  the  case  of  two  fasteners  at  an  angle  to  the 
load  direction,  a  nonconservative  failure  load  prediction  is  made. 

in  test  case  238,  a  70/20/10  laminate  layup  is  considered. 
The  fastener  loads  predicted  by  SAMCJ  correlate  well  with  the 
measured  values  in  Reference  2.  The  predicted  load  orientations  at 
fastener  locations  1  and  2  (see  Figure  41)  are  in  fair  agreement 
with  strain-gaged  bolt  measurements.  However,  the  measurement  at  tb 
third  location  is  suspect,  and  does  not  correlate  with  SAMCJ 
prediction.  SAMCJ  predicts  failure  at  11.4  kips,  in  a  shearout  mode 
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at  the  innermost  fastener  location.  Delaminations  accomompaniod  the 
predicted  shearout  failures  in  the  test  specimens  (Reference  2). 

The  predicted  failure  load  (11.4  kips)  is  12t  lower  than  the  average 
measured  value  (12.9  kips). 

SAMCJ  predictions  for  the  30/60/10  laminate  (test  case 
239)  are  presented  in  Figure  42.  The  fastener  loads  predicted 
correlate  well  with  the  measured  values  in  Reference  2.  The 
predicted  load  orientations  at  fastener  locations  1  and  2  are  in 
fair  agreements.  Again,  the  measurement  at  the  third  location  is 
suspect,  and  does  not  correlate  with  SAMCJ  prediction.  SAMCJ 
predicted  failure  at  11.7  kips,  with  a  net  section  failure  mode  at 
fastener  location  3.  Net  section  failures  occurred  at  both 
locations  number  3  and  2  in  the  test  specimens  (Reference  2) .  The 
predicted  failure  load  (11.7  kips)  is  lit  lower  than  the  average 
measured  value  (13.2  kips). 

3.5  Composite-to-Metal  Joints  with  Six  Fasteners  and  an 

Adjacent  Circular  Cut-Out  in  the  Laminate 

Test  cases  243,  246  and  247  address  a  single-shear  load 
transfer  between  0. 5-inch-thick  aluminum  plates  (without  a  cutout) 
and  40-ply  laminates  of  50/40/1.0,  70/20/10  and  25/60/15  layups,  with 
a  one  inch  diameter  circular  cut-out  adjacent  to  the  fasteners.  The 
bolted  laminates  with  the  circular  cutout  were  modeled  as  shown  in 
Figure  43,  for  the  three  test  cases. 

SAMCJ  predictions  for  the  40-ply  50/40/10  laminate  (test 
case  243)  are  presented  in  Figure  44.  SAMCJ  prediots  the  applied 
load  to  be  divided  nearly  equally  among  the  six  fasteners,  with 
fasteners  at  locations  2  and  5  carrying  the  largest  fraction.  This 
correlates  fairly  well  with  the  strain-gaged  bolt  measurements  in 
Reference  2.  SAMCJ  predicts  failure  to  occur  in  a  net  section  mode, 
across  the  1- inch-diameter  circular  cutout.  Two  out  of  three  test 
replicates  failed  in  the  predicted  manner.  One  replicate,  however, 
failed  in  a  net  section  mode,  across  the  inner  fasteners  (4,5,6). 


80 


E 


E 


3£ 


>2  Pj 


8AMCJ 

PREDICTION 

TE8T  RE8ULT8 
(Ref.  2) 

P1  /p 

0.26 

0.27 

P2/P 

0.41 

0.30 

P3/P 

0.33 

0.34 

<91  (degrees) 

14.0 

0.2 

0 2  (degrees) 

0.8 

-3.8 

Oq  (degrees) 

-12.0 

2.4 

Pf«llur.  (k">,> 

11.7 

13.2 

FAILURE  LOCATION 

3 

3,  2 

FAILURE  MODE(S) 

NET  SECTION 

NET  SECTION 

j 


P inure  43.  Nlnn-Element  Model  of  the  Bolted  Plates  in  Test  Cases  243, 
246  and  247 
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Figure  44.  SAMCJ  Predictions  and  Test  Results  for  Test  Case  243. 
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The  failure  load  predicted  by  SAMCJ  (38.3  kips)  is  9%  lower  than  the 
measured  average  value  (42.0  kips).  SAMCJ  predicts  a  higker  load 
level  (53.7  kips)  for  a  net  section  failure  across  the  in  ar 
fastener  holes ,  observed  in  one  out  of  three  test  replicaies. 

SAMCJ  predictions  for  the  40-ply,  70/20/10  lami1  ate  (test 
case  246)  are  presented  in  Figure  45.  Predicted  load  distribution 
among  the  six  fasteners  is  similar  to  that  predicted  for  the 
50/40/10  laminate,  and  is  in  fair  agreement  with  test  measurements. 
SAMCJ  predicts  a  shearout  mode  of  failure  at  an  inner  fastener 
location(5),  while  tests  resulted  in  shearout  failures  at  all  the 
fastener  locations  (1  to  6),  accompanied  by  delaminations.  The 
failure  load  predicted  by  SAMCJ  (37.9  kips  is  12%  lower  than  the 
measured  average  value  (42.9  kips).  It  is  also  noted  that  SAMCJ 
predicts  a  net  section  failure  across  the  1- inch-diameter  circular 
cut-out  at  a  slightly  higher  load  level  (38.8  kips). 

SAMCJ  predictions  for  the  40-ply,  25/60/15  laminate  with  a 
1-inch-diameter  circular  cut-out  (test  case  247)  are  presented  in 
Figure  46.  A  more  even  load  distribution  among  the  six  fasteners  is 
predicted  in  this  case,  and  is  in  agreement  with  strain- gaged  bolt 
measurements.  SAMCJ  predicts  failure  to  occur  at  30.6  kips,  in  a 
net  section  mode  across  the  circular  cutout  (7) .  Two  out  of  three 
test  replicates  in  Reference  2  failed  in  the  predicted  mode. 

However,  the  third  replicate  failed  in  a  net  section  mode  across  the 
inner  fastener  holes  (4,5,6).  SAMCJ  predicts  this  to  be  the  second 
probable  failure  at  a  higher  load  level  (36.0  kips).  The  predicted 
failure  load  (30.6  kips)  is  18%  lower  than  the  measured  average 
value  (37.2  kips). 

3.6  Composite-to-Metal  Joints  with  Five  Fasteners  in  a  Row 

Test  cases  250  and  251  consider  load  transfer  between  an 
aluminum  plate  and  a  40-ply,  50/40/10  laminate,  in  double-  and  single- 
shear  configurations,  respectively.  Figure  47  presents  the  six- 
element  model  of  the  bolted  plates  used  by  SAMCJ  for  both  test 
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Figure  46.  SAMCJ  Predictions  and  Test  Results  for  Test  Case  247. 


cases . 


SAMCJ  predictions  for  teat  case  250  are  presented  in 
Figure  48.  The  predicted  load  distribution  among  the  fasteners 
qualitatively  follows  the  trend  indicated  by  the  strain-gaged  bolt 
measurements  in  Reference  2.  The  predicted  peak  fastener  load 
fraction  (0.293),  however,  is  lower  than  the  measured  value  (0.345). 
The  innermost  fastener  (5)  is  predicted  to  carry  this  peak 
fractional  load,  in  agreement  with  test  results.  SAMCJ  predicts  the 
laminate  to  fail  in  a  net  section  mode  across  the  innermost  fastener 
hole  (5) .  This  prediction  is  also  in  agreement  with  the 
observations  in  Reference  2.  The  analytically  predicted  failure 
load  (13.4  kips)  is  24%  lower  than  the  measured  average  value  (17.7 
kips) . 

In  a  single-shear  load  transfer  situation  (test  caso  251) , 
a  more  even  load  distribution  among  the  five  fasteners  is  predicted 
by  SAMCJ  (see  Figure  49) .  This  agrees  well  with  the  straingaged 
bolt  measurements  in  Reference  2.  The  analytically  predicted  and 
experimentally  observed  failure  is  a  net  section  failure  across  the 
innermost  fastener  hole  (5) .  The  failure  load  predicted  by  SAMCJ 
(12.4  kips)  is  25%  lower  than  the  measured  average  value  (16.6 
kips) . 
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Test  Case  251,  Static  Tension,  Single-Shear 
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SECTION  4 


CONCLUSIONS 

A  strength  analysis  was  developed  for  laminates  that  are 
bolted  to  other  laminates  or  metallio  plates  by  many  fasteners.  The 
analysis  was  programmed  to  be  the  SAMCJ  computer  code.  The  validity 
of  the  developed  analysis  was  established  by  considering  different 
composite-to>metal  joints  that  were  tested  in  Reference  2.  SAMCJ 
accurately  predicted  the  observed  failure  location  and  the  failure 
mode  for  all  the  test  cases.  Predicted  fastener  load  distributions 
were  in  agreement  with  strain-gaged  bolt  measurements.  Predicted 
failure  loads  were  in  reasonable  agreement  with,  and  lower  than,  the 
measured  average  values  in  most  of  the  considered  test  cases. 
Nonoonservative  failure  loads  (larger  than  the  measured  value)  were 
predicted  only  in  a  few  test  cases.  In  summary,  the  developed 
strength  analysis  adequately  predicted  the  failure  load,  the  failure 
location  and  the  failure  mode  for  bolted  laminates. 

The  primary  limitations  of  the  developed  analysis  include 
its  inability  to  aacount  for  the  effect  of  countersunk  fasteners, 
its  inability  to  predict  the  precipitation  of  delaminations,  and  the 
inaccuraces  introduced  by  the  five-node  representation  of  a  complex 
problem  (a  laminate  with  finite  planform  dimensions,  a  fastener 
hole,  and  a  fastener  load  distribution  around  the  hole  boundary) . 

The  fastener  analysis  segment  of  SAMCJ  can  bo  modified  to  overcome 
the  first  limitation.  The  failure  procedure  can  be  modified  to 
predict  delaminations  through  approximate  estimations  of 
interlaminar  stresses  and  an  appropriate  failure  criterion.  This 
task  will  be  similar  to  that  performed  in  Reference  1.  The  last 
limitation  (inaccuracies  introduced  by  the  five-node  element)  can  be 
overcome  by  developing  a  higher  order  element  (nine  node  element) 
following  the  procedure  described  in  Section  2.  The  suggested 
improvements  of  the  developed  analysis  were  beyond  the  scope  of  this 
program,  and  are  recommended  as  future  efforts. 
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Despite  the  above  limitations,  the  developed  analysis 
(SAMCj  code)  offers  the  user  the  following  advantages:  (l)  SAMCJ  is 
a  test -independent,  one-step  analysis  that  oomputes  the  load 
distribution  (magnitude  and  orientation}  among  many  fasteners,  and 
subsequently  predicts  the  failure  load  as  the  lowest  value  of  the 
computed  applied  load  levels  for  each  of  three  failure  modes  at 
every  fastener  and  cutout  location;  (a)  SAMCJ  accounts  for  stress 
concentration  interaction  among  fasteners  and  adjacent  cutouts;  (3) 
SAMCJ  can  approximately  account  for  the  effect  of  tapering  of  bolted 
plates;  and  (4)  SAMCJ  is  a  validated  analytical  tool  that  oan  be 
used  in  the  dosign  of  efficient  joints  in  laminated  structural 
parts . 


Zn  summary,  the  strangth  analysis  of  bolted  laminates  (the 
SAMCJ  computer  code) ,  developed  in  this  Northrop/AFWAL  program,  is  a 
signifioant  contribution  to  tha  design  and  analysis  of  bolted 
laminated  structural  parts.  Its  current  limitations  do  not  restrict 
its  immediate  applicability  to  design  situations.  The  analysis  in 
moderately  conservative  and  provides  the  user  with  a  fairly  accurate 
prediction  of  the  failure  location  and  the  overall  failure  mode.  As 
such,  it  will  be  very  useful  in  rapidly  and  analytically  evaluating 
many  bolted  joint  concept a  ,  to  eeleot  the  most  efficient  concept 
for  a  defined  application. 
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G.  DOBEN 

Chief,  Scientific  &  Tech  Info  Gp 
Information  Services  Branch 


cc:  AFWAL/FIBRA 
(V.  Venkayya) 


UNITED  STATES  AIR  FORCE 


SEPTEMBER  18,1947 


